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We study the extremely polarized two-component Fermi gas with a mass imbalance in the strongly 
interacting regime. Specifically we focus on the experimentally available mixture of ^Li and ''"K 
atoms. In this regime spin polarons, i.e., dressed minority atoms, form. We consider the spectral 
function for the minority atoms, from which the lifetime and the effective mass of the spin polaron 
can be determined. Moreover, we predict the radio-frequency (rf) spectrum and the momentum 
distribution for the spin polarons for experiments with ®Li and ^''K atoms. Subsequently we study 
the relcixation of the motion of the Fermi polaron due to spin drag. 
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I. INTRODUCTION 

In many condensed-matter systems the response to a 
single impurity determines the low-temperature behav- 
ior of the system. Probably the most famous example 
hereof is a single electron moving in a lattice. Local lat- 
tice distortions, i.e., phonons, interact with the electron 
and together they form a quasiparticle that is known as 
the polaron because of the local change in polarization 
[l|. Another well-known impurity problem is that of an 
immobile magnetic impurity in a metal causing an en- 
hanced resistance below a certain temperature due to the 
Kondo effect 0]. The multichannel version of this effect 
has especially received much interest in the pastjbecause 
it leads to the formation of a non- Fermi liquid [3[. 

Here we study an impurity problem in a two- 
component atomic Fermi gas. An important motivation 
to use ultracold atoms is the unprecedented experimental 
control in these systems. They offer the interesting possi- 
bility of not only changing, for instance, particle numbers 
and temperature, but also the interaction strength. Via 
a Feshbach resonance the bare interaction can be tuned 
all the way from being weakly attractive (BCS-regime) 
to strongly attractive (BEC-regime), where in the inter- 
mediate regime the scattering length is much larger than 
the average interparticle distance. This so-called unitar- 
ity or strongly interacting limit, is the regime we focus 
on in this paper. 

We consider a mixture at zero temperature consist- 
ing of two (spin)species of fermions, where there is one 
minority particle immersed in a non-interacting sea of 
majority particles. The mass-balanced Fermi gas with 
high spin polarization has been studied extensively, both 
experimentally 0-@ and theoretically 0-[i3| ■ At the uni- 
tarity limit the minority particle gets dressed by a cloud 
of majority particles forming a quasiparticle similar to 
the polaron. This quasiparticle is often referred to as a 
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spin polaron, because its formation is due to interactions 
between particles in different spin states, or as a Fermi 
polaron, because it consists of fermionic atoms. Recently, 
the imbalanced spin-dipole mode the radio-frequency 
(rf) spectrum of the spin polaron [J|, and its energy and 
effective mass ^ that are different from those of the bare 
minority particle, have all been measured in this case. 

An intriguing new possibility for experiments is hav- 
ing a mass imbalance between the minority and the ma- 
jority particles by mixing two different atom species. 
A very promising mixture in this respect is the one of 
^Li and ""^K atoms. These atoms together have already 
been trapped and cooled to quantum degeneracy [15| . 
and moreover, several Feshbach resonances were identi- 
fied [l^. Theoretically, the phase diagram of the ^Li- 
''"K mixture has been determined [13, [3: and it dif- 
fers greatly from the phase diagram of a spin-imbalanced 
mixture by having not only a superfluid but also a super- 
solid region, depending on the sign of the polarization. 
We show here that already the two limiting cases of this 
mixture, i.e., a single light impurity in a sea of heavy 
atoms and vice versa, turn out to differ qualitatively in 
a manner that reflects the underlying asymmetry of the 
phase diagram. 

Indeed, in the solely spin-imbalanced case, having a 
cr) or a I — (t) minority particle results in the same im- 
purity problem, while with two different atom species 
there are two fundamentally different impurity problems. 
Thus, by introducing a mass imbalance, not only does 
the question whether dressed impurities still represent 
the ground state of the system arise, but so does the 
question of what is the difference between a heavy and a 
light impurity. Here, because the different atom species 
act as a pseudospin, the same many-body mechanism 
causes the dressing of the minority atom as for the solely 
spin-imbalanced case. Therefore, we also call this quasi- 
particle a spin polaron. In this paper we study for the 
two mass-imbalanced cases both a molecular bound state 
and the spin polaron. We show that although it does not 
form the ground state, the molecular bound state virtu- 
ally plays an important role in the system. In addition. 
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we study dissipation of kinetic energy of the minority 
particle due to interactions with the majority cloud that 
lead to spin drag [l9l |. 

In this paper we consider a homogeneous gas of atoms, 
while experiments are always done in a trap. Still, when 
1/kp is much smaller than the size of the cloud, where 
kp is the Fermi momentum of the majority atoms, the 
gas can locally be considered homogeneous and all our 
results apply. In this manner the appropriate averaging 
over the trap can be fully taken into account. 



II. MOLECULAR BOUND STATE 
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In the unitarity limit the minority particle interacts 
strongly with the Fermi sea of majority particles. Due to 
the low densities in ultracold atomic systems two-body 
processes represent the dominant scattering mechanism, 
where the minority particle can scatter off a majority 
particle an arbitrary number of times. Taking this into 
account in diagrammatic language results in an infinite 
sum of ladder diagrams, the so-called ladder sum. For 
the extremely imbalanced case at unitarity the bare in- 
teraction with the complete ladder sum added, i.e., the 
many-body T matrix, obeys 
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with Tvp and fiD. denoting the total momentum and en- 
ergy of the two incoming particles, and t {\) denot- 
ing a majority (minority) particle. The distribution 
for the majority particles is the Fermi-Dirac function 
N{x) = \l\e^l^°^ 1], with T the temperature and ks 
Boltzmann's constant. The distribution for a single mi- 
nority atom can be taken equal to zero in the integrand. 
The energy of an atom in state \a) is ^o-.k = Eo-.k — Mo-, 
where e^-.k = fi'^k'^ /2m„ and /ig- are the kinetic en- 
ergy and the chemical potential. The kinetic energy 
£k — (et.k + £i,k)/2 is associated with twice the reduced 
mass. Throughout this paper we take equal to the 
Fermi energy ep, while the chemical potential of the mi- 
nority atom is determined self-consistently from the self- 
energy [H, = KT.^{0, 0), defined later in Eq.©. 

A pole in the T matrix corresponds to a bound state, 
where the real part of the location of the pole is its en- 
ergy and the imaginary part is inversely proportional to 
its lifetime. In the above many-body T matrix the pole 
physically corresponds to a Feshbach molecule dressed 
by the Fermi sea, which we here refer to as a molecu- 
lar bound state. The energy Em{p) of this bound state 
at zero temperature, divided by the majority particles 
Fermi level ep, is shown as a solid line in Fig. [1] Up 
to some momentum Pmax it is a stable molecular state, 
while for larger momenta the imaginary part is non-zero 
and thus the bound state acquires a finite lifetime. A mi- 
nority and a majority atom cannot scatter off each other 
if their combined energy lies below a certain level due to 



FIG. 1: Dispersions of the molecular bound states e = huj+^i 
at r = 0, scaled by the Fermi energy ef ~ (Gn^n-f)'^^^ /2mf 
of the majority cloud, as a function of the momentum hp = 
ft|p| of the molecule, scaled by hkF ~ ^'ZmfEp- The grey area 
in both panels is the continuum of particle-particle excitations 
above the Fermi sea. The solid lines are the dispersions of the 
stable (decaying) molecular state when it lies below (in) the 
continuum. Panel (a) corresponds to one ®Li atom in a sea of 
^'^K atoms, panel (b) to one '*''K atom in a Fermi sea of ®Li 
atoms. 



Pauli blocking of the Fermi sea. Above this energy level 
there is a continuum of scattering states. This continuum 
of particle-particle excitations is also depicted in Fig. [TJ 
From the molecular dispersions it already becomes 
clear that a light impurity is very different from a heavy 
one. Fitting the dispersion of the molecular state for 
small momentum to Em{p) = ^'^P^ /'^'m-M^'rEMiO), shows 
that for the light impurity the stable molecule has a neg- 
ative effective mass, ttim — — 0.13m^, and has an energy 
Em{0) — ^Asp. The dispersion is qualitatively the same 
as for the mass-balanced case, where the stable molecule 
also has a negative mass, namely vriM — — 3.9m^. In 
contrast to the light impurity, with a heavy impurity 
the stable molecular state has a positive effective mass 
niM — 0.96m^ and has an energy Em{0) — —0.2ep. In- 
terestingly, it is also the part of the phase diagram with a 
minority of heavy particles that differs qualitatively from 
the mass-balanced case and contains a supersolid phase 
[TtI . [isj . In all cases the continuum of particle-particle 
excitations pushes the molecular state down, which is 
a consequence of level repulsion as in the more simple 
case of an avoided crossing of two energy levels. For the 
light impurity this repulsion results in a negative effective 
mass for the molecule. For the heavy impurity the effec- 
tive mass is positive, but smaller than one would obtain 
in the absence of the continuum. 



III. SPIN POLARON 

The presence of a molecular bound state does not mean 
necessarily that a molecule is the ground state of the sys- 
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FIG. 2: Spectral functions pi{e) of the minority particles at 
zero temperature and with zero momentum as a function of 
the energy e = hco + ^i- Panel (a) depicts the light impurity, 
panel (b) the heavy impurity. The delta peak in both figures 
corresponds to a stable polaron with energy Ep. 

tem, because some other state can have a lower energy 
than the molecule. We therefore now consider the dressed 
impurity, the spin polaron, and compare its energy with 
the molecule to determine the ground state of the system. 
The energy and lifetime of the quasiparticle can be ob- 
tained from the spectral function pi = — Im[G4,]/7r, where 
G_i is the Green's function describing the minority parti- 
cle in the presence of the Fermi sea. To obtain the latter 
a self-energy is added to the bare inverse Green's func- 
tion via G^^ — Gq I ~ ^i- At zero temperature and in 
the many-body T matrix or ladder approximation, that 
has been very successful for the mass-balanced case [lo| . 
we have 

n^^{q,u;+)^ y __Tk+q^^,^++^^ _^7V(e^,k), (2) 

with cli+ = lo + iO. Because the relevant momentum of 
the minority particle at zero temperature is much smaller 
than the Fermi momentum of the Fermi sea, we take its 
momentum equal to zero first. Then the spectral func- 
tion, for both impurity problems, has at the energy Ep 
a delta- function peak (see Fig. [2|), which corresponds to 
the energy of a stable quasi-particle, i.e., the spin po- 
laron. After comparing this energy with the energy of 
the molecular state Em{0) we conclude that for both 
cases the quasi-particle has lower energy and thus forms 
the ground state of the system. 

Apart from the energy of the dressed particle, also the 
quasi-particle residue Zp and the effective mass m* can 
be determined from the spectral function. The quasi- 
particle residue is the weight of the delta peak and the 
effective mass can be obtained from the momentum de- 
pendence of its location. For the light polaron, a dressed 
^Li atom, we find Ep ~ —2.2eF, Zp ~ 0.8 and m* ~ 
1.25mi, while for the dressed atom Ep ~ —OAAep, 
Zp ~ 0.64 and m* ~ l.lGm^. The energies and effective 
masses are in good agreement with previous theoretical 
results and Monte-Carlo calculations 10, 20] that do not 
consider the full spectral function. 

The presence of the molecular pole is very important 
for the spectral functions pj^(k, a;). In particular, the 
threshold of the continuum of pj^(0,w) is at zero energy 
when the molecular state always has a positive energy, as 



for the light impurity, see Fig. [2ja). In contrast, for the 
heavy impurity the molecular state can have a negative 
energy and this causes the threshold of the continuum 
to lie at a negative energy, see Fig. [DJb). The spectral 
function at zero temperature can be approximated by 
Piik,uj) ~ ZpS{el + Ep-fii-huj), with = h^k^/2m*. 
For both impurity problems, however, it does not capture 
all the features of p|(k, w), as we will see next. 

A direct probe for the quasiparticle residue Zp is the 
momentum distribution of the minority particles, which 
can be obtained experimentally by a time-of-flight exper- 
iment. From the spectral function it can be calculated 
by means of A'^;(fc) = / dujpi{k,uj)N{oj). In Fig. [3] the re- 
sults are shown, both for the full spectral function and for 
the delta peak only, at zero temperature and for polariza- 
tion P — 0.9. Also depicted are the ideal gas momentum 
distributions for comparison. The quasiparticle residue 
can be read off easily in both figures. It can also be 
seen that the delta peak is a good approximation for the 
heavy impurities, while for the light impurities Zp de- 
pends more strongly on the external momentum, which 
is not captured by this approximation. 

The energy of the spin polaron can be directly ob- 
tained from the rf spectrum, which was recently mea- 
sured for the mass-balanced case. In an rf experiment 
incoming photons with frequency w^/ induce transitions 
from an occupied hyperfine state to an empty state. 
The fraction of transferred atoms as a function of the 
photon frequency is the rf spectrum, where the thresh- 
old of the spectrum is the polaron energy. Theoreti- 
cally, the spectrum can be obtained directly from the 
spectral function by using the Kubo formula, I{uJrf) oc 
/ dk7V(^;,k - hujrf)Pii'k,^lM - ^rf) m. When using 
the low-temperature spectral function the integral can be 
performed, yielding 

IK,) oc ZpMiu:., + Ep)N f "^^"-f + ) . (3) 

V \ m* — mi J 

The rf spectra for the two mass-imbalanced impurity 
problems are shown in Fig. |4] for P = 0.99 and the 
temperature of the experiment with mass balance, T = 
0.1 ^4]. For the light impurity the analytic result 
from Eq. ^ reproduces the full spectral function result 
almost exactly. 

IV. SPIN-DRAG RELAXATION RATE 

At zero temperature the spin polaron corresponds to 
a delta-function peak in the spectral function as we have 
just seen. At nonzero temperatures we expect this peak 
to broaden and to obtain a width that is proportional 
to at low temperatures. An immediate consequence 
of this non-zero width is that the polaron acquires inter- 
esting transport properties. In particular it leads to a 
non-zero spin-drag relaxation rate 1 /Tgd of the polaron 
moving in a Fermi sea of majority particles. The friction 
of the spin polaron and the out-of-phase dipole mode 
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are examples of properties determined by l/r^d. For the 
mass-balanced case the latter has been studied experi- 
mentally [5, Gij, and the transport properties of the mass- 
imbalanced impurity problem have been studied theoret- 
ically using thermodynamic arguments to calculate the 
effective interaction [19l |. 

We derive here a general expression for the relaxation 
rate of one polaron moving with velocity v through a 
cloud of majority particles with which it interacts, where 
its velocity is small compared to the Fermi velocity of 
the majority particles |v| <C kp/m^Ti. The equation of 
motion of the spin polaron then reads 
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FIG. 4: RF spectra at T = Q.IATf and polarization P = 0.99 
in arbitrary units. The solid lines are obtained from the full 
spectral functions, the dashed line from Eq. ([Sj. In panel (a) 
the results for the light impurity are depicted and in panel 
(b) for the heavy impurity. 



^ ^ ^r(v) ^ 

at m*ni T^d 



(4) 



where r(v) is the Boltzmann collision integral, which was 
linearized in the last step. For the spin-drag relaxation 
rate for the impurity problem we obtain in this manner 
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where /? is I/fcsT and A'^(k) is the distribution function 
of the majority particles. The on-shell effective interac- 
tion Vic,k',q in general depends on the incoming momenta 
k and k' and on the transferred momentum q of the scat- 
tering particles. From the linearized collision integral, 
the above expression is obtained by using as the distri- 
bution function for the dressed impurity a delta function. 
The result in Eq. ([5]) is generic for any impurity, fermionic 
or bosonic, in any environment, fermionic or bosonic. 

In the impurity problem at hand we take for iV-|-(k) the 
Fermi-Dirac distribution function. At low temperatures 
only small q contribute and the difference between the 
two distributions becomes strongly peaked around the 
Fermi level [2l|. For the interaction we take the many- 
body T matrix from Eq. ([T]) , with an additional factor 
Zp to account for the wavefunction renormalization, and 



then ultimately obtain 
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Tsd 12TThm*^ (27r) 



tP-\T{kF,eF)\^ [ dq- 

T^r J SI 



sinh^(/3e*/2) 



(6) 



where "/{m^/m-^) is a dimensionless function depending 
on the mass ratio of the minority and majority parti- 
cles. For the light impurity we find 7(0.15) ~ 8.58, while 
we find 7(6.7) ~ 1.96 for the heavy impurity. The tem- 
perature dependence for only one minority particle in a 
fermionic environment is the same as for the spin-drag 
relaxation rate for equal densities of fermions, namely 
1/Tsd oc T^. This quadratic temperature dependence is 
expected for a Fermi liquid and recently it was verified 
experimentally for the mass-balanced case that l/rgd in- 
deed decreases as the temperature decreases [1, [Ij] • The 
result in Eq. ^ implies that at T = there is no spin 
drag relaxation, which in turn implies that the spin po- 
laron is a stable quasiparticle in that case. As mentioned 
above, the latter is consistent with the delta-peaks in 
the spectral functions in Fig. |4] and confirms that the 
ladder approximation captures the relevant physics for 
these mass-imbalanced mixtures. 
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FIG. 3: Momentum distributions N{k) for the spin polarons 
at zero temperature and P — 0.9, where the polarization P = 
{n-f — ni)/{nf + ni). The solid, black (gray) lines are the 
momentum distributions of the polarons obtained from the 
full (delta-peak) spectral function, while the dashed lines are 
the distributions for an ideal gas. Panel (a) depicts the light 
spin polaron and panel (b) the heavy one. 



V. CONCLUSION 

We calculated a number of important observables of 
the extremely polarized ^Li-'^'^K mixture. We showed 
that at the unitarity limit, although virtually the molec- 
ular state plays an important role, polarons form at low 
temperatures and dominate all physical properties of the 
mixture. Apart from its equilibrium properties we also 
looked at the transport properties of the spin polaron and 
found that the spin-drag relaxation rate takes a universal 
form and scales with the square of the temperature, as 
expected for a Fermi liquid. 
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